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Abstract 



We consider the swelling kinetics of layered structures. We focus on the case 
of triblock copolymer mesogels, although our results are applicable to other 
layered structures including clays. We assume the mesogels are swollen by 
a solvent that is good for the bridging block but poor for the non-bridging 
block. At long times the penetration front moves as in ordinary diffusion, i.e., 
as t 1 / 2 . At short times, however, the bending elasticity of the non-bridging 
layers becomes important. This bending elasticity leads to a i 1 / 6 relaxation 
of the penetration front at early times. The crossover length between these 
two regimes is approximately the width of a single layer. However, for a 
large number of lamellae there is a cooperative effect which leads to a large 
enhancement of this crossover length. 
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I. INTRODUCTION 



In the strong segregation limit, a symmetric diblock copolymer with A and B blocks will 
form a lamellar phase consisting of a series of AB layers. This kind of system has received 
much study and may have applications as a novel composite material. In general, however, 
it suffers from one major drawback: the AB lamellae are only weakly bonded to one another 
by local Van der Waals forces, and possibly by a few entanglements. The system thus lacks 
mechanical integrity in the melt state and the AB interfaces are very weak. This is the case 
even if both blocks are glassy. If for instance the system is placed in a solvent which is good 
for the B blocks, these blocks swell and the individual lamellae separate from one another. 
However, as noted by Halperin and Zhulina [1], one very simple modification allows much 
greater mechanical integrity. This is to use an ABA triblock copolymer. A certain fraction 
of the B chains bridge between layers. In the melt state this system, though a liquid, is 
much more difficult to pull apart, because of the presence of these bridges. If the A regions 
are rubbery or glassy then the system becomes a solid, even though the B regions may still 
show liquid-like behavior. When the B regions are swollen by a selective solvent (which 
is good for B and bad for A), the system becomes a sandwich with rubbery liquid layers 
between solid regions. This swollen system was given the name "mesogel" [1]. Although 
the theoretical studies of these systems have been rather recent, mesogels have existed in 
practice for much longer [2] . Their mechanical integrity offers many advantages over diblock 
systems, and they are already being used in controlled drug release systems [3]. 

Previous theoretical studies have concentrated on the equilibrium behaviour and to a 
lesser extent on the rheology of mesogels [1]. The questions of interest for equilibrium 
properties are the degree of bridging, the degree of swelling, and the mechanical moduli of 
the composite system. Here our interest lies in the kinetics of swelling, i.e., the approach 
to equilibrium. The swelling kinetics of isotropic gels is a subject for which there is a 
broad literature reporting many novel effects [4-6]. Mesogels, are by definition non-isotropic, 
and might exhibit interesting swelling behaviour because of this as well as the competition 
between layer elasticity and solvent penetration. The swelling of lamellar mesogels is related 
to the swelling of other layered systems such as clays [7]. These have numerous applications, 
particularly in environmental science. Geometrically the two systems are similar, although 
there are major differences in both the type and magnitude of the physical forces involved. 
In clays, for instance, the layers are much smaller and electrostatic and short-range forces are 
thought to be of great importance. Another related area is the swelling of polymeric glasses. 
In that case there are two kinds of temporal behaviour, ordinary diffusion and "diffusion" 
linear in time. There have been many proposals for why this occurs [8,9], all of which agree 
more-or-less with the experimental results. However, at present there is no widely accepted 
definitive theory, although elastic effects almost certainly play some role. 

In any system involving mesogels the gels must at some time go through a swelling 
process, and our study is applicable to all such cases. The opposite case where one de-swells 
a mesogel may be important in drug delivery systems [10]. The same model can be used to 
describe de-swelling, although we do not discuss it further here. In section II we examine the 
case of a single layer of a mesogel, such as might be obtained by spin-coating. This turns 
out to have a mesoscopic length L±, below which the elasticity of the layers needs to be 
accounted for, and for which the swelling kinetics are non-diffusive in character. It turns out 
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that L\ is of order the layer thickness. In section III we generalise this to more macroscopic 
multilamellar samples. There, the crossover length becomes much larger: Lm — M 1 ^ 2 Li, 
where M is the number of layers. Hence Lm can become of macroscopic size. We conclude 
with a discussion in section IV. 

II. THE MODEL OF DIFFUSION PROCESS IN SINGLE LAYER 

Consider a melt of ABA triblocks, with blocks of roughly equal length. At low temper- 
atures this system will microphase separate to form a series of distinct B layers separated 
by A layers, with sharp boundaries in between. In general this will have a large number 
of defects, which can be removed by shearing. At lower temperatures, or if the A phase is 
crosslinked, the A region may become a rubbery solid. The system then consists of a series 
of rubbery A layers separated by molten B layers. This results in a well-ordered "meso- 
rubber". This system can be placed in a selective solvent that is good for the B block and 
poor for the A block. We assume for simplicity that only one edge of the specimen is in 
contact with the solvent reservoir. The solvent penetrates the B regions and swells them. 
Here we consider a single ABA lamella. Initial contact of the solvent with the B regions 
is energeticly favourable but is opposed by two effects. The first is the stretching of the 
B chains to form a brush. The second is the bending and stretching of the surrounding A 
regions. This implies that the free energy per unit area of the system consists of three terms 



Here we ignore the stretching energy of A domains, which is asecond-order effect in relaxed 
layers. We need to evaluate each of these terms in the partially swollen state to describe the 
swelling kinetics. 

Before swelling, the B regions form a dense melt of thickness h. After equilibrium swelling 
is reached, the B regions form a swollen polymer brush [11,12] with the same grafting 
density as the melt. In the swollen state, dense grafting enforces strong overlap among 
the undeformed coils. For tethered chains in a good solvent, this increases the number 
of monomer-monomer contacts and the corresponding interaction energy. This penalty is 
reduced by stretching the chains along the normal to the grafting sites, thereby lowering the 
monomer concentration in the layer and increasing the layer thickness h. Stretching lowers 
the interaction energy per chain, Fb interaction! at the price of a higher elastic free energy, 
Fb stretch- The interplay of these two terms sets the equilibrium thickness of the layer h eq . 

We write the layer thickness of the B region in the partially swollen state as 



where ip is the solvent concentration in the B region. We can express the volume fraction of 
solvent (f) s and the volume fraction of polymer <p p as <p s = + ip) and (f) p — 1/(1 + ■0). In 
the Alexander-de Gennes model, the B layer is envisioned as a close-packed array of blobs 
of uniform size £, which is equal to the grafting distance d. Each of these blobs costs a free 
energy of kT. The free energy of the B layer per unit area is given by [11] 



F Fb stretch ~\~ Fb interaction ~\~ Fj± bend- 



(1) 



h(x) = h(l + ip{x)), 



(2) 
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Fb Fb interaction Fb stretch 



kT[^)^+kT(-J^-]4, (3) 



<r>) '" ' " \Na 2 d 2 ) 



where numerical prefactors have been ignored, and a is the B monomer size. Since <p p 
(1 + ip)' 1 and h = h(l + if)), Eq. (|3]) can be rewritten as 



h 



a 3 



(JL)n 



F B = kT | - ] (1 + i-)-' + kT [^j^ j (1 + i-)-- I ' ) 



Minimising this yields the equilibrium swelling ip eq = (hNd 2 /lafij 3 — 1. It is convenient to 

simplify Fb by expanding about ip eq , to yield an approximate free energy per unit area for 
the B regions [13,14] 

F fl ~ - ^) a . (5) 

Here, the coefficient S in Eq. (||) is given by 

o d 2 F B 
dtp 2 

[Na 2 d 2 ) [~NcP 

The remaining term in the free energy is the bending elastic energy of the A regions: 

F A =\n{v"h) 2 = \k' (7) 

where k is the bending elastic constant and k' = h 2 K. If the A regions form a rubbery 

solid, then we can estimate k from elasticity theory as k ~ n r H\ where fi r is the modulus 
of rigidity of the A regions and is the thickness of the A layers [15]. This modulus is 
related to the number of crosslinks per unit volume n, by fi r ~ nkT. 
The chemical potential of the solvent is 

5F 

= K'(v 4 ij)+S(iP-ij eq ). (8) 

The solvent current is defined by 3 S = — A^V/i, where A4 is a mobility, which we take to 
be independent of the concentration ip [16]. From the continuity equation 

^ + V-J s = 0, (9) 
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we have the diffusion equation 



f = v.,^) 

= k'M (v 6 </>) + SM (V 2 ^) . (10) 

Here we consider a single layer parallel to the xy plane, with infinite extent in the y direction 
(Fig. H). The solvent lies in the region x < and the gel in the region x > 0. This reduces 
the problem to one dimension, with spatial gradients only in the x direction. 

From Eq.(|H]), we expect the diffusion process to be dominated by the bending energy 
at shorter lengthscales (or earlier times), and to be dominated by the interaction energy at 
longer lengthscales (or later times). Therefore, the diffusion process is characterized by a 
crossover time scale T\ and a crossover length L\. Depending on the initial conditions, for 
times less than T\ or for solvent penetration less than L 1; the process is not characterized 
by simple diffusion. Instead, the profile spreads as t 1 / 6 . At later times, as the solvent profile 
spreads out, simple diffusion dominates. In most of what follows, we shall be concerned 
primarily with the non-diffusive regime, and in particular with crossover length, below which 
it can be observed. We can estimate the crossover length L\ by the substitution V — > \jL\. 
From Eq. ([H]), the crossover length is approximately 

In a rubbery solid where k ~ kTnH\, the crossover length L\ becomes 



nH 3 A Na 2 d 2 



a . 



3 



(12) 



Within the Alexander-de Gennes approximation, the chains are assumed to be equally 
stretched. The grafting distance is estimated to be 
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\ r y a2 J 

We can also estimate the layer thickness of A domains as 

Na 3 



d~ (1^ | Nt>a. (13) 



H 



A 



AkT 



i 

' 3 ,2 



Nza. (14) 



The layer thickness of the B domains, prior to swelling is similarly given by h ~ aN 2 ^ 3 . 
Substitution of Eqs. (|13|) and (JTJ) into Eq. Ql2D leads to a crossover length 
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Lx ~ (na 3 y/ 4 N 13 / 72 h 



(15) 



If we consider a tightly crosslinked network so that na 3 ~ 1, then the crossover length is 

roughly the layer thickness h. This is a mesoscopic length of order 10 2 to 10 3 Angstroms. 
Note, however, that if the modulus of the A regions is very large (say, for instance, if the A 
regions were glassy), then the crossover length L\ could be much larger. Before the solvent 
has reached x = Li, we expect the solvent front moves as t 1//6 . At later times it obeys 
ordinary diffusion — moving as t 1 ^ 2 . The anomalous t 1//6 "diffusion" is a result of the bending 
elasticity of the A layers. A similar anomalous exponent can be found in the hydrodynamics 
of membrane systems [17,18]. 



III. THE MODEL OF DIFFUSION PROCESS IN MULTILAYERS 

In the above we considered the swelling of a single layer. In that case, the A layer 
elasticity only affected the initial swelling up to a length scale of order of the layer thickness. 
For a multilayer system, this effect can be much more pronounced. As illustrated in Fig. |2|, 
the swelling of the inner layers is transmitted to the outer layers, and hence the outer layers 
must bend significantly more. In order to examine this, we first consider the case of uniform 
swelling, in which each layer is swollen by the same amount at a given value of x. Below, we 
shall relax this assumption. For M layers, uniform swelling corresponds to ip m (x) = ip{x) 
for m = 1, 2, • • • M. The height of m-th layer is given by 

h m (x) =mh(l + ip(x)) . (16) 
The free energy of the mesogel is 

M r l /__,. M 1_ . 



m=l 



M 3 , / 2 ,N2 SM 



_ K '(vV) + — ^-^ eq y. (17) 



Thus 



?t ~ K ' M Ml V 6 iP + SMMV 2 ^. (18) 
at 3 

In other words, the relaxation rate for a mode of wavevector q in the x direction, is 

V 1 ~ M + SMq 2 ^ . (19) 

The characteristic length Lm, beyond which the swelling becomes dominated by ordinary 

diffusion, is found by equating the two terms on the right hand side of (|19"D with q = 1/ Lm, 
and is 
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L m kM*L x , 



(20) 



where L\ is the crossover length for a single layer. For a large number, M, of layers, this is 
a lengthscale much larger than the thickness of a single lamella. However, since we consider 
the case of swelling that occurs from an edge, our results are valid for a mesoscopic crossover 
length Lju less than the horizontal dimensions of the sample in Fig. 2. The corresponding 
crossover time is 

r M ~ Mh 2 /D, (21) 

where D is the diffusion constant for the solvent. 

To go beyond the approximation of uniform swelling, we express the height of the m-th 
layer in a more general form 

h m (x) - ft, m _! (x) = h(l-if) m (x)), (22) 

where m — 1, 2, • • • M, and h = 0. The free energy of the mesogel is 

1 M ( „ q 

z m=l k ' t 

The chemical potential of Z-th layer can be calculated as 



h m - h m -l - k(l + Ipeq) 



(23) 



M 

E 



\ Sh m (x) dipt (x) 

M g M-l 

= nllY^ V 4 /i m + j- ( 2h m ~ h m+1 - 

m=l m=l 

The dynamical equations for hi are given by 



h m -i) + 



h M - h M -i - h{l + ipeg) . (24) 



d(hi - hi-!) 
dt 



M 



M-l 



Mh 2 K V 6 /i m + MS V 2 (2hm - h m+1 - hm-x) 

n = 

+ MSV 2 (h M - h M -i) 



and 



d(h M - h M -i) 
dt 



nMh 2 V & h m + MSV 2 (h M - h M -i) 



(25) 



(26) 



For large M and small displacements of the layer heights away from their initial positions, 
we characterize the layer heights by a displacement field u (x, z): 



h m (x) = mh + u (x, mhj , 



(27) 



Eq. ( P5| ) becomes 
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9 ( dhl \ - -MkV'u Or, z) + MSV 2 ^ V 



dt \ dz 2 



dz 2 



(28) 



after differentiating once with respect to z. (A discrete expression of this can also be obtained 
by subtraction of Eq. |25j for adjacent /.) If we express the displacement field u(x, z) as a 
sum over wavevectors in the x and z directions, q and k z , 



u( Xl z) = Y.u(q } k z )e^ x+k * z \ 
q,k z 



then the relaxation rate function of q and k z is 

t-- 1 = MKq 6 k; 2 + MSq 2 



(29) 



(30) 



The initial swelling mode corresponds to the smallest value of k z consistent with the bound- 
ary conditions, namely 



k 7 



1 



Mh 



(31) 



Thus, the crossover length is 



i.e. the same result as was found in the uniform swelling case (p0|). 

Equation (p3|) can also be expressed in terms of the concentration variables ip m (x) 



(32) 



where 



l M r 2 

1 m=l ^ 



■ 9^ r • 



(33) 



= + ^ ( X )) 

1=1 

The chemical potential of l-th layer is given by 



(34) 



M m 

Hi = Kh 2 V "Vm' + S (ipl - Ipeq) 

m=l m'=l 

The dynamical equations for ipi can be expressed as 

9 -'• x / M ■ ■ ■ 2 1 \ / V 6 ^i \ 

Kh 2 M 



V J 



2 ■•■ 2 1 

\i ••• i ly 



V 6 Vm-i 

v v 6 ^Af y 



v v 2 ^Af y 



(35) 



(36) 



S 



At early times, for which the spatial gradients are large, the dominant relaxation rate cor- 
responding to a mode of wavevector q is 



t" 1 = Kh 2 MX {max) q 6 + SMq 2 , (37) 

where A ( - max < 1 is the largest eigenvalue of the matrix 

( M ■■■ 2 1\ 

1 . (38) 

2 ••• 2 1 

\1 •■• 1 l) 

It can be shown that the eigenvalues of this matrix are given by [19] 

1 

2 (1 — cos a; 



A = 7777^ (39) 



where uj satisfies 

t ;H ! ( \ f,.:' 1 f( it ^ 



tan(Mcu) = cot^. (40) 



Thus, for large M, the leading eigenvalue is 



A (max) ^ 4 m2 (41) 



7T 2 V M 

(This was also determined by numerical evaluation of A ( - max ) for M up to 100.) Thus, the 
crossover length 

L U * Ll M^ (l + JL) . (42) 

On simple physical grounds it is expected that the anomalous swelling at early times will 
be most apparent for the inner layers, since the curvature elasticity of the outer layers tends 
to constrain the inner layers. In other words, the initial relaxation of a profile such as the 
one shown in Fig. 2 will be dominated by a mode that relaxes gradients of the inner layers 
(/ii(x), h 2 (x), . . .) faster than the outer layers (. . . , h M _i(x), h M (x)). This is reflected in the 
eigenvector corresponding to this mode (Eq. pll), which is given by 



(max) . ( {M-l + l)u \ 
W j ocsm^ — j. (43) 

Thus, indeed the initial relaxation of the inner layers (small I) is greater than that of the 
outer layers (large I). However, the crossover length in Eq. ( |42"D still grows as M 1//2 , as for 
uniform swelling. 
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IV. DISCUSSION 



In this paper we have examined the swelling kinetics of lamellar mesogels. For such 
mesogels the elasticity of the solid lamellae leads to a novel kinetics which is non-diffusive 
at short times. Similar effects should be important for the swelling of clays and other 
layered structures. There is one better-known case of non-diffusive behaviour, and that 
occurs in the swelling of isotropic polymer glasses [8,9]. In that system there occurs both 
case I behaviour, which obeys ordinary diffusion, and case II behaviour where the front 
moves linearly in time. There is even a "super-case II" behaviour which has a front position 
moving as t 3//2 . In general the explanations of case II involve either a diffusion constant which 
strongly depends on concentration, or a stress relation which similarly depends strongly on 
the concentration. Lamellar mesogels represent a totally different type of behaviour, where 
the front moves as t 1 / 6 , at early times and as t 1 ! 2 at later times. Note that here we have 
ignored the dependence of the diffusion constant on concentration and have expanded the 
free energy about the equilibrium state. Our results are thus most applicable to the situation 
where the meso-rubber is first pre-swollen and then swollen again to reach an equilibrium 
swelling. However, the general effect of the A layer elasticity will always be to modify the 
swelling at early times. We have also ignored cracks in our analysis. If the strains set up in 
the system are very large then cracks may appear in the A region. Such cracks are indeed 
seen in the swelling of polymer glasses, which can swell explosively. 
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FIGURES 



FIG. 1. The swelling geometry for a single layer. The solvent penetrates from the left and the 
gel is assumed infinite towards the right and out of the plane of the page. The swelling profile is 
characterized by h(x). 

FIG. 2. The swelling geometry for multiple layers. Note that in this case the swelling of the 
inner layers causes a large bending at the outer layers. 
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